The RHESSI Microflare Height Distribution

We present the first in-depth statistical survey of flare source heights observed by RHESSI. Flares were found using a flare-finding algorithm designed to search the 6-10 keV count-rate when RHESSI's full sensitivity was available in order to find the smallest events (Christe et al., 2008). Between March 2002 and March 2007, a total of 25,006 events were found. Source locations were determined in the 4-10 keV, 10-15 keV, and 15-30 keV energy ranges for each event. In order to extract the height distribution from the observed projected source positions, a forward-fit model was developed with an assumed source height distribution where height is measured from the photosphere. We find that the best flare height distribution is given by g(h) \propto exp(-h/{\lambda}) where {\lambda} = 6.1\pm0.3 Mm is the scale height. A power-law height distribution with a negative power-law index, {\gamma} = 3.1 \pm 0.1 is also consistent with the data. Interpreted as thermal loop top sources, these heights are compared to loops generated by a potential field model (PFSS). The measured flare heights distribution are found to be much steeper than the potential field loop height distribution which may be a signature of the flare energization process

Nonlinear force-free reconstruction of the global solar magnetic field: methodology

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday's equation, give rise to a respective normal field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modelled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition - the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.Comment: 18 pages, 5 figures, Solar Physic

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. II. Numeric Code and Tests

Based on a second-order approximation of nonlinear force-free magnetic field solutions in terms of uniformly twisted field lines derived in Paper I, we develop here a numeric code that is capable to forward-fit such analytical solutions to arbitrary magnetogram (or vector magnetograph) data combined with (stereoscopically triangulated) coronal loop 3D coordinates. We test the code here by forward-fitting to six potential field and six nonpotential field cases simulated with our analytical model, as well as by forward-fitting to an exactly force-free solution of the Low and Lou (1990) model. The forward-fitting tests demonstrate: (i) a satisfactory convergence behavior (with typical misalignment angles of $\mu \approx 1^\circ-10^\circ$), (ii) relatively fast computation times (from seconds to a few minutes), and (iii) the high fidelity of retrieved force-free $\alpha$-parameters ($\alpha_{\rm fit}/\alpha_{\rm model} \approx 0.9-1.0$ for simulations and $\alpha_{\rm fit}/\alpha_{\rm model} \approx 0.7\pm0.3$ for the Low and Lou model). The salient feature of this numeric code is the relatively fast computation of a quasi-forcefree magnetic field, which closely matches the geometry of coronal loops in active regions, and complements the existing {\sl nonlinear force-free field (NLFFF)} codes based on photospheric magnetograms without coronal constraints.Comment: Solar PHysics, (in press), 25 pages, 11 figure

The Magnetic Field of the Solar Corona from Pulsar Observations

We present a novel experiment with the capacity to independently measure both the electron density and the magnetic field of the solar corona. We achieve this through measurement of the excess Faraday rotation due to propagation of the polarised emission from a number of pulsars through the magnetic field of the solar corona. This method yields independent measures of the integrated electron density, via dispersion of the pulsed signal and the magnetic field, via the amount of Faraday rotation. In principle this allows the determination of the integrated magnetic field through the solar corona along many lines of sight without any assumptions regarding the electron density distribution. We present a detection of an increase in the rotation measure of the pulsar J1801$-$2304 of approximately 160 \rad at an elongation of 0.95$^\circ$ from the centre of the solar disk. This corresponds to a lower limit of the magnetic field strength along this line of sight of $> 393\mu\mathrm{G}$. The lack of precision in the integrated electron density measurement restricts this result to a limit, but application of coronal plasma models can further constrain this to approximately 20mG, along a path passing 2.5 solar radii from the solar limb. Which is consistent with predictions obtained using extensions to the Source Surface models published by Wilcox Solar ObservatoryComment: 16 pages, 4 figures (1 colour): Submitted to Solar Physic

Multiwavelength Study on Solar and Interplanetary Origins of the Strongest Geomagnetic Storm of Solar Cycle 23

We study the solar sources of an intense geomagnetic storm of solar cycle 23 that occurred on 20 November 2003, based on ground- and space-based multiwavelength observations. The coronal mass ejections (CMEs) responsible for the above geomagnetic storm originated from the super-active region NOAA 10501. We investigate the H-alpha observations of the flare events made with a 15 cm solar tower telescope at ARIES, Nainital, India. The propagation characteristics of the CMEs have been derived from the three-dimensional images of the solar wind (i.e., density and speed) obtained from the interplanetary scintillation data, supplemented with other ground- and space-based measurements. The TRACE, SXI and H-alpha observations revealed two successive ejections (of speeds ~350 and ~100 km/s), originating from the same filament channel, which were associated with two high speed CMEs (~1223 and ~1660 km/s, respectively). These two ejections generated propagating fast shock waves (i.e., fast drifting type II radio bursts) in the corona. The interaction of these CMEs along the Sun-Earth line has led to the severity of the storm. According to our investigation, the interplanetary medium consisted of two merging magnetic clouds (MCs) that preserved their identity during their propagation. These magnetic clouds made the interplanetary magnetic field (IMF) southward for a long time, which reconnected with the geomagnetic field, resulting the super-storm (Dst_peak=-472 nT) on the Earth.Comment: 24 pages, 16 figures, Accepted for publication in Solar Physic

TomograPy: A Fast, Instrument-Independent, Solar Tomography Software

Solar tomography has progressed rapidly in recent years thanks to the development of robust algorithms and the availability of more powerful computers. It can today provide crucial insights in solving issues related to the line-of-sight integration present in the data of solar imagers and coronagraphs. However, there remain challenges such as the increase of the available volume of data, the handling of the temporal evolution of the observed structures, and the heterogeneity of the data in multi-spacecraft studies. We present a generic software package that can perform fast tomographic inversions that scales linearly with the number of measurements, linearly with the length of the reconstruction cube (and not the number of voxels) and linearly with the number of cores and can use data from different sources and with a variety of physical models: TomograPy (http://nbarbey.github.com/TomograPy/), an open-source software freely available on the Python Package Index. For performance, TomograPy uses a parallelized-projection algorithm. It relies on the World Coordinate System standard to manage various data sources. A variety of inversion algorithms are provided to perform the tomographic-map estimation. A test suite is provided along with the code to ensure software quality. Since it makes use of the Siddon algorithm it is restricted to rectangular parallelepiped voxels but the spherical geometry of the corona can be handled through proper use of priors. We describe the main features of the code and show three practical examples of multi-spacecraft tomographic inversions using STEREO/EUVI and STEREO/COR1 data. Static and smoothly varying temporal evolution models are presented.Comment: 21 pages, 6 figures, 5 table

Beyond the Singularity of the 2-D Charged Black Hole

Two dimensional charged black holes in string theory can be obtained as exact (SL(2,R)xU(1))/U(1) quotient CFTs. The geometry of the quotient is induced from that of the group, and in particular includes regions beyond the black hole singularities. Moreover, wavefunctions in such black holes are obtained from gauge invariant vertex operators in the SL(2,R) CFT, hence their behavior beyond the singularity is determined. When the black hole is charged we find that the wavefunctions are smooth at the singularities. Unlike the uncharged case, scattering waves prepared beyond the singularity are not fully reflected; part of the wave is transmitted through the singularity. Hence, the physics outside the horizon of a charged black hole is sensitive to conditions set behind the past singularity.Comment: 19 pages, 5 figures; v2: refs added, minor typos corrected; v3: references on the infinite blue shift at the inner horizon and minor corrections adde

Multiwavelength Study of M8.9/3B Solar Flare from AR NOAA 10960

We present a multi-wavelength analysis of a long duration white-light solar flare (M8.9/3B) event that occurred on 4 June 2007 from NOAA AR 10960. The flare was observed by several spaceborne instruments, namely SOHO/MDI, Hinode/SOT, TRACE and STEREO/SECCHI. The flare was initiated near a small, positive-polarity, satellite sunspot at the centre of the AR, surrounded by opposite-polarity field regions. MDI images of the AR show considerable amount of changes in a small positive-polarity sunspot of delta configuration during the flare event. SOT/G-band (4305 A) images of the sunspot also suggest the rapid evolution of the positive-polarity sunspot with highly twisted penumbral filaments before the flare event, which were oriented in the counterclockwise direction. It shows the change in orientation and also remarkable disappearance of twisted penumbral filaments (~35-40%) and enhancement in umbral area (~45-50%) during the decay phase of the flare. TRACE and SECCHI observations reveal the successive activations of two helical twisted structures associated with this sunspot, and the corresponding brightening in the chromosphere as observed by the time-sequence images of SOT/Ca II H line (3968 A). The secondary-helical twisted structure is found to be associated with the M8.9 flare event. The brightening starts 6-7 min prior to the flare maximum with the appearance of secondary helical-twisted structure. The flare intensity maximizes as this structure moves away from the AR. This twisted flux-tube associated with the flare triggering, is found to be failed in eruption. The location of the flare is found to coincide with the activation site of the helical twisted structures. We conclude that the activations of successive helical twists in the magnetic flux tubes/ropes plays a crucial role in the energy build-up process and triggering of M-class solar flare without a CME.Comment: 22 pages, 12 figures, Accepted for Publication in Solar Physic

Deflection and Rotation of CMEs from Active Region 11158

Between the 13 and 16 of February 2011 a series of coronal mass ejections (CMEs) erupted from multiple polarity inversion lines within active region 11158. For seven of these CMEs we use the Graduated Cylindrical Shell (GCS) flux rope model to determine the CME trajectory using both Solar Terrestrial Relations Observatory (STEREO) extreme ultraviolet (EUV) and coronagraph images. We then use the Forecasting a CME's Altered Trajectory (ForeCAT) model for nonradial CME dynamics driven by magnetic forces, to simulate the deflection and rotation of the seven CMEs. We find good agreement between the ForeCAT results and the reconstructed CME positions and orientations. The CME deflections range in magnitude between 10 degrees and 30 degrees. All CMEs deflect to the north but we find variations in the direction of the longitudinal deflection. The rotations range between 5\mydeg and 50\mydeg with both clockwise and counterclockwise rotations occurring. Three of the CMEs begin with initial positions within 2 degrees of one another. These three CMEs all deflect primarily northward, with some minor eastward deflection, and rotate counterclockwise. Their final positions and orientations, however, respectively differ by 20 degrees and 30 degrees. This variation in deflection and rotation results from differences in the CME expansion and radial propagation close to the Sun, as well as the CME mass. Ultimately, only one of these seven CMEs yielded discernible in situ signatures near Earth, despite the active region facing near Earth throughout the eruptions. We suggest that the differences in the deflection and rotation of the CMEs can explain whether each CME impacted or missed the Earth.Comment: 18 pages, 6 figures, accepted in Solar Physic

Automated Coronal Hole Detection using Local Intensity Thresholding Techniques

We identify coronal holes using a histogram-based intensity thresholding technique and compare their properties to fast solar wind streams at three different points in the heliosphere. The thresholding technique was tested on EUV and X-ray images obtained using instruments onboard STEREO, SOHO and Hinode. The full-disk images were transformed into Lambert equal-area projection maps and partitioned into a series of overlapping sub-images from which local histograms were extracted. The histograms were used to determine the threshold for the low intensity regions, which were then classified as coronal holes or filaments using magnetograms from the SOHO/MDI. For all three instruments, the local thresholding algorithm was found to successfully determine coronal hole boundaries in a consistent manner. Coronal hole properties extracted using the segmentation algorithm were then compared with in situ measurements of the solar wind at 1 AU from ACE and STEREO. Our results indicate that flux tubes rooted in coronal holes expand super-radially within 1 AU and that larger (smaller) coronal holes result in longer (shorter) duration high-speed solar wind streams